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ABSTRACT 

We apply the adhesion approximation to study the formation and evolution of 
voids in the Universe. Our simulations - carried out using 128 3 particles in a cubical 
box with side 128 Mpc - indicate that the void spectrum evolves with time and that 
the mean void size in the standard COBE-normalised Cold Dark Matter (hereafter 
CDM) model with h$o = 1, scales approximately as D(z) = r^ppg ) where Do ~ 10.5 
Mpc. Interestingly, we find a strong correlation between the sizes of voids and the 
value of the primordial gravitational potential at void centers. This observation 



could in principle, pave the way towards reconstructing the form of the primordial 
potential from a knowledge of the observed void spectrum. Studying the void spec- 
trum at different cosmological epochs, for spectra with a built in /c-space cutoff we 
find that, the number of voids in a representative volume evolves with time. The 
mean number of voids first increases until a maximum value is reached (indicating 
that the formation of cellular structure is complete), and then begins to decrease as 
clumps and filaments merge leading to hierarchical clustering and the subsequent 
elimination of small voids. The cosmological epoch characterizing the completion of 
cellular structure occurs when the length scale going nonlinear approaches the mean 
distance between peaks of the gravitational potential. A central result of this paper 
is that voids can be populated by substructure such as mini-sheets and filaments, 
which run through voids. The number of such mini-pancakes which pass through 
a given void, can be measured by the genus characteristic of an individual void 
which is an indicator of the topology of a given void in initial (Lagrangian) space. 
Large voids have on an average a larger genus measure than smaller voids indicat- 
ing more substructure within larger voids relative to smaller ones. We find that 
the topology of individual voids is strongly epoch dependent, with void topologies 
generally simplifying with time. This means that as voids grow older they become 
progressively more empty and have less substructure within them. We evaluate the 
genus measure both for individual voids as well as for the entire ensemble of voids 
predicted by the CDM model. As a result we find that the topology of voids when 
taken together with the void spectrum, is a very useful statistical indicator of the 
evolution of the structure of the Universe on large scales. 
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I. INTRODUCTION 



A remarkable feature of the large scale structure of the Universe is the presence 
of voids - volumes of order 10 3 h~ 3 Mpc 3 , virtually devoid of the presence of galaxies. 
The existence of voids was spectacularly demonstrated with the discovery of the 
Bootes void an ~ 10 5 /i~ 3 Mpc 3 underdense region, by Kirshner et al. (1981). Since 
then successive redshift surveys (de Lapparant et al. 1986, Vogeley et al. 1991; 
Slezak et al. 1993), as well as deep pencil beam surveys (Broadhurst et al. 1990), 
have confirmed that voids are a salient feature of the large scale structure of the 
Universe, and that galaxies seem to lie preferentially along sheets and filaments 
separating voids. 

Although it is generally agreed that voids account for most of the volume of 
the Universe, there is still considerable disagreement in the literature as to what is 
the typical size of a void. Part of this disagreement has to do with how one chooses 
to define voids. If voids are defined as underdense regions then the size of the void 
will obviously depend upon the density threshold below which an underdense region 
"becomes" a void. Clearly, the mean volume of underdense regions will always be 
larger than the mean volume of completely empty regions. The average size of 
a void in a sample will, in addition, be sensitive to the smallest voids which one 
recognizes. Clearly, if we recognize empty regions of diameter ~ 0.2 Mpc to be 
voids, then the mean void size will be much smaller than if we were dealing with 
voids having diameters > 2 Mpc. 
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However, even if we restrict ourselves to voids defined as completely empty re- 
gions, we still come across differing points of view as to the typical void size in the 
Universe. Thus, Kauffmann and Fairall (1991) (henceforth KF), who consider only 
void diameters larger than ~ 5/i _1 Mpc in an uncontrolled data sample, find from 
their catalogue of over 100 voids that the typical (viz, modal) size of a completely 
empty region is D ~ 10/i -1 Mpc. The mean void diameter in the KF sample is 
36.5/i _1 Mpc (Little 1992). 1 The authors also find no evidence of a cutoff in void 
size upto a maximum void diameter of 64 Mpc. On the other hand, Vogeley et 
al. (1991), who work with a smaller but controlled sample (the CfA survey), find 
the maximum size of a completely empty region in their sample to be ~ 20h~ 1 Mpc 
which is considerably smaller than the mean void size in the KF sample. A part of 
this discrepency may arise from the fact that KF work with a large but uncontrolled 
data sample, whereas the CfA survey although containing fewer galaxies has the 
advantage of being magnitude limited. (KF use galaxy redshifts which are not mag- 
nitude limited since they are compiled from many different sources, thereby giving 
rise to the "distinct possibility that certain regions of the sky are undersampled in 
comparison to other regions " and to the related possibility that the existence of a 
particular void is due to undersampling rather than to any physical process (Kauff- 
mann & Fairall 1991).) The other reason why the results of Vogeley et al. (1991) 
and those of KF differ, could be due to the different techniques used to actually 

1 Little has evaluated the mean void diameter for 129 voids determined from 
the merged Southern Redshift Catalogue and the Catalogue of Radial Velocities of 
Galaxies and listed by KF in their paper. 
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measure void sizes. Thus, Vogeley et al. (1991) construct a void probability func- 
tion for the CfA survey which measures the probability that a randomly selected 
volume in their survey contain no galaxies, whereas KF determine a void spectrum 
by filling regions empty of galaxies in their sample with cubes, the edge length of 
the largest cube filling a given empty region providing a measure of the void size. 
Although this method leads to a well defined algorithm to measure void sizes it 
suffers from the drawback (acknowledged by the authors) of ignoring the issue of 
void topology. (Voids in the KF sample are assumed a-priori, to possess a bubble 
like topology.) 

Since voids are associated with positive peaks in the primordial gravitational 
potential, (Gurbatov, Saichev & Shandarin, 1985, 1989; Kofman & Shandarin, 
1988) they give rise to anisotropies in the Cosmic Microwave Background Radi- 
ation (CMBR) via the Sachs-Wolfe effect: (ST/T) ~ 5(p/3 on the surface of last 
scattering. This effect has been studied for individual voids by Blumenthal et al. 
(1992) who find the sizes of the largest voids to be smaller than ~ 130/i _1 Mpc (one 
void of this size within the horizon). In a later paper the same authors incorporate 
the COBE results to give a stronger limit ~ 60/i _1 to the maximal void size (Piran 
et al. 1993). 

From the theoretical viewpoint the fact that voids might provide a key to un- 
derstanding the large scale structure of the Universe was emphasised more than a 
decade ago, by Zeldovich & Shandarin (1982) and Zeldovich, Einasto & Shandarin 
(1982). Since then, the evolution of voids has been examined both analytically 
as well as numerically by a number of authors (Peebles 1982; Hausman, Olson 
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& Roth 1983; Hoffman, Salpeter & Wasserman 1983; Fillmore & Goldreich 1984; 
Bertschinger 1985; Blaes, Villumsen & Goldreich 1990; Bonnor & Chamorro 1990; 
Regos & Geller 1991; Dubinski et al. 1992; van de Weygaert & van Kampen 1993; 
Harrington, Melott & Shandarin 1993). Semi-analytical spherically symmetric stud- 
ies of voids have shown that large voids can often arise out of small initial negative 
density fluctuations. Numerical simulations by Bertschinger 1985 and Blaes, Vil- 
lumsen & Goldreich 1990, confirmed that aspherical negative density perturbations 
grow to become more spherical with time (Icke, 1984) (the time-reverse of the Lin, 
Mestel & Shu (1965) instability), indicating that the evolution of negative density 
perturbations can often be treated as the time reversal of positive density ones. 
Since \Sp/p\ never exceeds unity within a void, the linear approximation might be 
expected to hold for a longer period inside a void. Also, as voids tend to expand 
at faster rates than the mean Hubble flow of the Universe, matter within voids will 
have a tendency to to be swept up into two dimensional sheets separating neighbor- 
ing voids. Both these considerations have led to the construction of a geometrical 
model of large scale structure based on Voronoi tesselation, in which voids act as 
centers of repulsion and matter collects in sheets, filaments and nodes, which to- 
gether make up the skeleton of the large scale structure of the Universe (Icke & van 
de Weygaert 1987; van de Weygaert 1991). 

In this paper we propose to study the properties of voids using a semi-analytic 
approach to model non-linear gravitational instability, known as the adhesion ap- 
proximation (Gurbatov, Saichev & Shandarin 1985, 1989; Shandarin & Zeldovich 
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1989). Since the adhesion model incorporates the Zeldovich approximation, we ex- 
pect the predictions of the latter to be also valid for the former at early times. Most 
voids, according to the Zeldovich approximation, can be successfuly associated with 
peaks in the primordial gravitational potential. However, at late times, after shell 
crossing and the consequent breakdown of the Zeldovich approximation, the inher- 
ently nonlocal nature of the adhesion model begins to manifest itself. At such times, 
the location and size of a void is determined not so much by the location and height 
of peaks, as by the global structure of the primordial gravitational potential. As a 
result the association of peaks in the gravitational potential with voids, which exists 
in the Zeldovich approximation, begins to progressively breakdown at late times. 
The adhesion model, like the model based on Voronoi tesselation, describes the 
Universe in terms of a skeletal structure consisting of pancakes (sheets) , filaments, 
clumps and voids. However, barring this resemblance, important differences exist 
between the two models. Namely, voids are not present abinitio in the adhesion 
model but arise out of dynamical considerations. Furthermore, the rich dynamical 
structure of the adhesion model allows it to address a number of important issues 
concerning voids, such as: Do voids evolve with time ? Are the voids seen today 
primary (i.e., first generation) or were they formed through the merger of an earlier 
generation of voids ? What is the likelyhood that voids will be populated by sub- 
structure such as mini-pancakes ? What is the topology of the large scale structure 
of the Universe ? etc. We shall attempt to address some of these questions plus a 
few more in the present paper. 

The organisation of our paper is as follows: 

7 



We briefly discuss the adhesion approximation in section 2, and compare some 
of its predictions with those of the Zeldovich approximation, and N-body simula- 
tions. In section 3 we apply the adhesion model to obtain a distribution of void 
sizes - the void spectrum, for idealised primordial spectra P(k) oc /c n , n = — 1, 0, +1 
(Melott & Shandarin 1993), as well as for more realistic initial conditions attempt- 
ing to explain the large scale structure of the Universe, such as the cold dark matter 
model. We find that the void spectrum evolves with time - the mean void diameter 
increasing as the Universe expands. We also find that a distinct correlation exists 
between the mean void size and the value of the primordial gravitational potential 
at the void center - larger voids being associated with higher regions of the pri- 
mordial potential. Our results also show that, in several instances a void can be 
multiply connected indicating a nontrivial topology, with minor pancakes (or fila- 
ments) running through a void that is circumscribed by major Zeldovich pancakes. 
The void topology evolves with time, proceeding from a sponge like topology at 
early times, to a bubble like topology at late times. Our results concerning the void 
topology are summarised in Section 4. A discussion of our results is presented in 
section 5. 

II. THE ADHESION MODEL AND THE ZELDOVICH APPROXIMATION 

The adhesion model is a logical extension of the Zeldovich approximation, which 
describes the evolution of density perturbations in a collisionless, self-gravitating 
medium, until the epoch of caustic formation and shell crossing (Zeldovich 1970; 
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for review see Shandarin & Zeldovich 1989). We assume a spatially flat, matter 
dominated Universe: fi m = fi = 1, A = 0. 

Particle trajectories in the Zeldovich approximation, are described by the map- 
ping 

x(q,t) =q-a(t)V$ (1) 

where q is the initial (i.e., Lagrangian) comoving coordinate, and x is the final (i.e., 
Eulerian) comoving coordinate of a particle. a(t) is the scale factor which coin- 
cides with the growing mode of the gravitational instability in a matter dominated 
spatially flat Universe. $ is the linear velocity potential related to the primordial 
Newtonian gravitational potential (f> by 

( p = A- 1 <S> (2) 

where A = 2/(3H 2 a 3 ) 7 H is the Hubble parameter: H 2 oc a~ 3 implying A is a 
constant. The Poisson equation 

^ = Aa(t)V 2 (t> (3) 
P 

together with the linear growth law dp/p oc a(t), implies that the linearised New- 
tonian potential 0, like its counterpart the velocity potential $, does not depend 
upon time. 

The continuity equation p(x)d s x = p(q)d s q (p(x) is the Eulerian density, and 
p(q) = P = constant is the Lagrangian density), when combined with the Zeldovich 
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approximation (1) describes the deformation of an infinitesimal volume element in 
terms of the Jacobian of the transformation from q to x, so that 



dV L 

and, as a result 



dx 




dq 





5ij - a(t) 



dqidqj 



(4a) 



dV E = dV L [1 - a(t)X 1 (q)] [1 - a(t)X 2 (q)} [1 - a(t)X 3 (q)} (46) 

where 6Ve and dVi,, are the volume elements in Eulerian space (E-space) and 
Lagrangian space (L-space) respectively, and Ai, A2 and A3 are the eigenvalues of 
the deformation tensor ip^ : 



dqidq k 



From (4b) we see that one (or more) of the eigenvalues Aj have to be positive, in 
order for (4) to describe a contraction along one (or more) axis (i.e., dVs < dV£)- 
On the other hand all three eigenvalues Ai,A2,A3 must be negative, in order for 
the volume element to expand along all three axis (i.e., 6Ve > dVi,)- It is well 
known that in the Zeldovich approximation, the condition (dVE < dVi)-, ultimately 
gives rise to the formation of caustics such as pancakes, filaments and clumps, 
whereas the opposite condition (dVE > dVi) determines voids. The validity of 
(4) breaks down soon after the formation of caustics, which occur when a(t) = 
min{l/Xi}. Thereafter the thickness of pancakes grows unbounded in the Zeldovich 
approximation, which runs counter to the findings of N-body simulations which 
show that the thickness of pancakes remains considerably smaller than both the 
pancake size as well as the mean distance between pancakes. In order to stabilise 
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the thickness of pancakes, a mock viscosity term, which mimicks the effects of 
nonlinear gravity, must be incorporated into the Zeldovich approximation. 

To see how this may be done, consider the Euler equation 

dv 1 1 

— + -(vV)v + Hv = — (5) 
at a a 

where 

dx 

v = a irv 

x being the comoving coordinate of the particle (or fluid element) . In terms of the 
new comoving velocity variable 

dx v 
da aa 

(5) can be recast as (Kofman 1991) 

^ + { uV)u = -^-{u + AV<t ) ) (6) 
aa la 

The right hand side of equation (6) is reminiscent of a force term, in the Lagrangian 
approach. Setting it to zero we obtain the Zeldovich approximation. In the adhesion 
approximation, we replace the right hand side of equation (6) by a viscosity term 
which mimicks the adhesive effects of nonlinear gravity. As a result the Euler 
equation for the velocity field becomes 

^ + {uV)u = vV 2 u (7) 
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which is the Burgers equation (Burgers 1974). For potential flows equation (7) has 
the following exact analytical solution 

f d 3 q ^ x (S) exp (J^G(q, x, a)) 

U(x , a) = - a{t) K2v ' (8) 

Jd 3 q exp (±G(q,x,a)) 

where 

G(S2,a) = »(4)-£^, 

$>(q) is the potential of the linear velocity field introduced earlier: u{x,a = 0) = 
— V$. For small values of the viscosity parameter v — > 0, the integrals in equation 
(8) can be evaluated using the method of steepest descents. This leads to the 
following condition being imposed on the gradient of the gravitational potential: 

V*P(g; x , a) = V^(g) = AV g >(# (9) 

where 

Equation (9) has an elegant geometrical interpretation which can be summarised as 
follows: Whether or not a particle originally located at qo is stuck within a pancake 
can be found by descending a paraboloid P(q; x, a) with radius of curvature 2a(t) 
and height Po onto the linear gravitational potential 0, in such a manner so as to 
be tangential to the potential at qo (see Fig. 1). (The height of the paraboloid Po 
is a free parameter uniquely determined by the condition that P(q, x, t) osculates 
4>(q) at qo. We use cj> and $ interchangeably, since </> oc $ if the Universe is flat 
and matter dominated) . If the paraboloid so constructed touches or intersects the 
potential at any other point q ^ qo, then we say that the particle in question 
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has already entered a caustic (Fig. 1, middle and bottom panels), otherwise it has 
not (Fig.l, top panel). Thus, one can divide the Lagrangian space at any given 
time into stuck and free regions. Stuck regions correspond to particles already in 
caustics, whereas free regions are still expanding via the Zeldovich relations (1,4) 
and correspond to voids. The apex of the paraboloid x, which osculates the potential 
in two or more points (without intersecting elsewhere) describes the location of the 
caustic in Eulerian space. The mass in caustics (clumps) in Fig. 1 is given by 
m i = \qi ~ 92 1 wi2 = |?3 — ?4| Pi and M = \q[ — q' A \ p, p being the mean density. 
(The clumps are drawn so that the clump radius is proportional to its mass.) The 
velocities of the caustics can be found from u = dx/da = (cf) qi — 4> qj )/(qi — qj) x n, 
where the unit normal n points from a higher value of the gravitational potential 
towards a lower value. (For earlier implementations of the adhesion model using 
the geometric approach, see Kofman, Pogosyan & Shandarin 1990, Sahni 1991, 
Williams et al. 1991, Kofman et al. 1992. For an alternative implementation of the 
adhesion approximation see Nusser & Dekel 1990, and Weinberg & Gunn 1990.) 

The skeleton of the large scale structure can be derived from the division of 
L-space into stuck and free regions, by noting that the border between these regions 
consists of particles which are only just entering into caustics. Since the Zeldovich 
approximation is valid right until a particle ends up in a caustic, we can determine 
the precise location of caustics in E-space by moving the border between stuck and 
free Lagrangian regions by means of the Zeldovich approximation. We do this for 
different values of a{t) defined in units of the rms linear density contrast on the 
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grid scale: 




oc 



2 



a(t) = ((Sp/pf) 



P(k)k 



(10) 



where N D is the dimensionality of the space. Our results are shown in Fig. 2 
a-e. In Fig. 2a, b the distribution of caustics (dots) was obtained by moving the 
border between stuck (unshaded) and free (shaded) Lagrangian regions by means 
of the Zeldovich approximation (equation (1)). The area of clumps in Fig. 2c is 
drawn proportional to the clump mass obtained using the adhesion approximation. 
Comparison of the adhesion model (Fig. 2b), with the two-dimensional N-body 
simulations of Melott and Shandarin 1989 (Fig. 2d), shows very good agreement at 
an epoch when the Zeldovich approximation (Fig. 2e) breaks down. 

An interesting question concerns the percolation of stuck and free regions in 
L-space. From Fig. 2 a-c we clearly see that for small values of a (a < a\) the 
free phase percolates (Fig. 2a), whereas for large values of cr, (a > a 2)1 percolation 
occurs in the stuck phase (Fig. 2c) . It would be interesting to estimate the values 
of cri,o"2, and their dependence on the power spectrum P(k). It would also be 
interesting to determine whether an epoch exists, when both phases percolate (i.e., 
whether ct 2 < cri). These questions might shed some light on the issue of the 
large scale topology of the Universe (i.e., whether "bubble-like" or "sponge-like"), 
and will be examined in detail in a forthcoming paper (Sahni, Sathyaprakash & 
Shandarin 1993). 

The geometrical picture of the adhesion model, allows us to make some gen- 
eral statements regarding the relative concentration of matter in clumps, filaments, 
pancakes and voids, which are the structural units constituting the skeleton of the 
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large scale structure of the Universe in three dimensions. The whole process of 
structure formation can be viewed as a mapping of the space at the initial time 
when the density is virtually homogeneous (the Lagrangian space) into the space 
at the time in question (the Eulerian space) . All structural units except voids have 
at least one dimension e (thickness) which is much smaller than the other dimen- 
sions £ « I . e is determined by the relaxation processes and in the adhesion 
model is parametrised by the viscosity coefficient v. The limiting case v — > is 
accompanied by e — > (for the relation between e and v see Gurbatov, Saichev & 
Shandarin 1989) and, as a result, clumps become points, pancakes become sheets of 
infinitesimal thickness etc. In this case, at late times, the mapping becomes formally 
degenerate: clumps correspond to the mapping of three-dimensional volumes onto 
zero-dimensional points (3D — > OD), filaments correspond to the mapping of two- 
dimensional surfaces onto one dimensional lines (2D — > ID), pancakes correspond 
to the mapping of one-dimensional lines onto two-dimensional surfaces (ID — > 2D), 
and voids (OD -> 3D). 

The relative concentration of matter within these structural units can be qual- 
itatively characterized by two numbers Vl,Ve, specifying the type of mapping 
T>l — > T>e '■ T>e is the dimension of the object in E-space, which is formed from 
matter originally spread uniformly in a T>l dimensional volume in L-space. Of 
course, the mapping of a higher dimensionsional set onto a lower dimensional set 
(clumps and filaments) formally results in infinite density, however of different types. 
Clumps, corresponding to the 3D — > OD mapping are much higher concentrations 
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of mass than filaments which result from the 2D — > ID mapping. In addition, fila- 
ments originating from a 2D set in 3D space have infinitesimal mass which means 
that they eventually vanish in real systems. Voids and pancakes result from the 
OD — > 3D and ID — > 2D mappings respectively and therefore tend to have both 
infinitesimal densities and masses at later times. 

This discussion may be useful in providing a guideline for a qualitative under- 
standing of the evolution of the structural units with time as well as for comparing 
3-D and 2D systems. In two-dimensional space clumps correspond to the 2D — > OD 
mapping, filaments to the ID — > ID mapping and voids to the OD — > 2D mapping. 
Therefore one may expect filaments in two-dimensional simulations (ID ID) 
to appear relatively more conspicuous when compared to pancakes (ID — > 2D) in 
three-dimensional simulations. On the other hand filaments in three-dimensional 
simulations (2D — > ID) must look more noticeable than their counterparts in two 
dimensions. We believe that this type of reasoning is worth keeping in mind when 
extrapolating (even qualitatively) conclusions from 2D to 3D. 

III. THE VOID SPECTRUM AND ITS EVOLUTION 

It is worth contrasting the dynamical picture of voids which emerges in the ad- 
hesion model, with the essentially static picture that follows from the Zeldovich ap- 
proximation. In the Zeldovich approximation, voids are associated with Lagrangian 
regions in which all the eigenvalues of the deformation tensor are negative. This 
implies essentially, that there is a one to one correspondence between voids and 
peaks in the primordial gravitational potential. This is also true at early times in 
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the adhesion model when the paraboloid is very narrow and can only feel small 
localised regions of </> (see Fig. 1, top panel). The adhesion model at such times 
gives essentially the same results as the Zeldovich approximation. At late times 
however, the radius of curvature of the paraboloid can become comparable to that 
of the initial gravitational potential, with the result that the paraboloid is now able 
to survey larger regions of the potential. It is at this stage that the essentially 
nonlocal nature of the adhesion model emerges. As a result, a free region in L-space 
(a void) can become stuck at late times, even if it lies in the vicinity of a peak in 
the gravitational potential (see Fig. 1, bottom panel). As a result, no strictly one 
to one relationship exists between voids and peaks in the gravitational potential, in 
the adhesion model. In addition, voids in L-space tend to shrink in the adhesion 
approximation, whereas in the Zeldovich approximation, their volume in L-space re- 
mains constant. (By contrast the total volume of voids in E-space remains constant 
in the adhesion model.) The distinction between voids in the Zeldovich approxima- 
tion and in the adhesion model in two-dimensions, is brought out clearly in Fig. 3 
in which we superimpose voids in the Zeldovich approximation with those obtained 
using the adhesion model. From Fig. 3 we see that virtually all of the voids in 
the adhesion model, constitute a subclass of voids in the Zeldovich approximation, 
although this is strictly true only at very late times, when most of the matter is al- 
ready in pancakes, so that the volume (equivalently, area in 2D) occupied by voids 
in L-space, is very small. The terms early and late can be quantified in relation 
to the rms linear density contrast on the grid scale a(t) defined in equation (10), 
which evolves in proportion to the expansion factor of the Universe. (The adhesion 
picture plotted in Fig. 3, corresponds to a = 8.) 
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The geometrical ansatz which allows us to divide L-space into stuck and free 
regions, also provides us with a very elegant way to determine both the individual 
void size, as well as the entire void spectrum. As discussed in the previous para- 
graph, a void in L-space is a region, matter from which has not yet fallen into a 
caustic. Since the Zeldovich approximation is still valid in such regions, one can 
use the volume deformation formula given by equation (4b), to determine the Euler 
volume occupied by a unit cubical element in L-space. Summing over all such ele- 
mentary volume elements within a given Lagrangian void gives the volume of the 
void in E-space: 

N 

V E = Y J dV E = Y J dV L [1 - a^Aite)] [1 - a(t)X 2 (q t )] [1 - a(t)X 3 ( qi )} (11) 
i=i 

where N denotes the number of elementary free volume elements in a given L-space 
void. 

(Note that voids in L-space are defined on a discrete grid, whereas those in 
E-space are not. Consequently Vl = NoIVl = N(\n/2) 3 gives the volume of the 
void in L-space, where is the Nyquist wavelength.) 

This method, which is a useful synthesis of both the adhesion as well as the 
Zeldovich approximations, allows us to determine the precise volume of a void with- 
out resorting to any simplifying assumptions regarding its shape, such as sphericity, 
etc., though it neglects the volume occupied by clumps, filaments and pancakes. 
We apply this ansatz to determine the spectrum of void sizes, for a Cold Dark 
Matter model with Q = 1, H = 50 km/sec/Mpc, and amplitude normalised by 
the COBE observed CMBR anisotropy (The COBE - normalised amplitude gives 
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(AM/M) ~ 1.1 for the r.m.s. CDM mass fluctuation on scales of 8h 1 Mpc (Smoot 
et al. 1992, Efstathiou, Bond & White 1992). The corresponding value of a(t) is 
9.) 2 . 

The results of our simulations, which average three realisations of the CDM 
spectrum using 128 3 particles, in a 128 Mpc box with periodic boundary conditions, 
are shown in Fig. 4 for different cosmic epochs corresponding to: (a) the past 
a = 0.5 ao, z = 1 (top panel); (b) the present a = ao, z = (middle panel); and 
(c) the future (!) a = 1.5 ao (bottom panel); (a - is the expansion factor, and z the 
redshift) . We plot both the number fraction of voids (right panels) having a given 
diameter: n(D)/N, N = Xli n (-^i)? as weu as the associated void volume fraction 
(left panels): v(D)/V = n(D) x V E (D)/V, V = EiVeCA)- (The void diameter 
in our simulations is defined so that 4tt(D/2) 3 /3 = Ve{D), where Ve is the void 
volume in Eulerian space, described by equation (11).) The void spectrum in the 
middle panel of Fig. 4 has been plotted for ~ 3400 voids, - the mean number of voids 
at the present epoch in our simulation. We find that the volume spectrum of voids 
shows a fairly uniform spread for voids having diameters in the range 10 < D < 30 
Mpc where D is the void diameter. Approximately 65% of the total volume in voids, 
is contained in voids lying within this range. By contrast, the number spectrum 
of void sizes, has a well defined peak on scales of ~ 6 Mpc, indicating that the 

2 In normalising the CDM spectrum, we have ignored the gravity wave contri- 
bution to (AT/T), which in some cases can be significant (Davis et al. 1992, Liddle 
&Lyth 1992, Lidsey k Coles 1992, Lucchin et al. 1992, Sahni & Souradeep 1992, 
Salopek 1992, Souradeep & Sahni 1992). 
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most abundant voids in our simulation (~ 50%) are of this size. However, the net 
contribution from such small size voids, to the overall void volume is a tiny ~ 10%. 

From Fig. 4 we see a distinct evolution of the void spectrum with cosmological 
epoch, the average size of a void growing roughly with time as D(z) = + z, 

where D ~ 10.5 Mpc, is the mean diameter of a void today and z is the cosmological 
redshift. 3 (The modal diameter of a void in our simulations turns out be slightly 
larger: D mode ~ 17 Mpc. The maximum void size in our simulations is ~ 60 Mpc.) 

The void spectrum obtained by us using the adhesion model, may be compared 
with the void spectrum reconstructed by KF from a preliminary catalogue of 129 
voids . (The voids in our simulations are empty regions, the density of matter in 
a void being less than 1% of the mean density of matter in the Universe. Likewise 
the voids in the KF sample are also assumed to be empty.) 

KF obtain D mode ~ 8 — ll/i -1 Mpc, h = if/100, for the typical (viz modal) 
diameter of a void in their catalogue. Applying the same criteria as KF, and con- 
sidering only voids having diameters greater than 5/i _1 Mpc in our sample, we find 
the mean, median and modal void sizes in our simulations to be D mean = 19.1 ±0.3 
Mpc, Dmedian = 26.4 ± 0.8 Mpc, and D mo d e = 25 ± 5 Mpc respectively, at the 

3 The mean void diameter is evaluated using 4n(D/2) 3 /3 = V, where V = 
iV _1 Xlili i s the mean void volume. Taking an alternate definition for D : 
D = iV -1 Eili A results in a somewhat smaller value D ~ 7 Mpc. 
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present epoch. 4 As a result we find that the modal diameter of a void in our sample 
is just about the same as that in the KF sample. (However, the mean diameter 
of a void in our simulations appears to be smaller than that implied in the KF 
survey.) The good agreement shown between the observed modal void size and that 
predicted by the CDM model must however be treated with some caution, since (a) 
the modal void size in our CDM simulations is subject to fairly large uncertainties 
(in fact, we do not have a well defined modal value at all) and (b) the methods 
we employ to determine the void size and those employed by KF are somewhat 
different. Namely, KF work with an uncontrolled data set and may therefore be 
systematically overestimating void sizes. KF further make the assumption that 
voids are essentially compact so that two neighboring voids do not percolate, they 
then determine the void size by inscribing the given void with a cube - the edge 
length of the cube providing a measure of the void size. While this method has the 
advantage of providing a well defined implementational algorithm to determine void 
sizes, it suffers from the drawbacks that: (i) voids are assumed to be approximately 
ellipsoidal with comparable axes, and (ii) because voids are assumed to be isolated, 
the possible sponge-like topology of voids is ignored. Neither of these assumptions is 
made in the adhesion model, in which the void volume is calculated for voids which 

4 The mean, median and mode are defined with respect to the void volume 
fraction. Voids having the modal void size therefore provide the largest contribution 
to the overall void volume spectrum. Evaluating statistical averages for the void 
number spectrum, gives much smaller values both for the median as well as the 
modal void size due to the abundance of small voids in the CDM simulation. 
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can have virtually any shape and topology - a point that we shall elaborate on in 
the next section. Consequently, in order to rule out (or confirm) the predictions 
made by the CDM model, we must apply the same void finding algorithm both to 
the data as well as to the simulations. We propose to do this in a follow up paper. 

An important related issue which we have been unable to adequately settle 
in the present paper concerns the distribution of void sizes in a realistic (i.e., spa- 
tially unbounded) CDM Universe. Simulations performed by us on 64 3 Mpc 3 and 
128 3 Mpc 3 sample volumes of a CDM Universe show that, although larger voids 
are present in the 128 3 sample relative to the 64 3 one, the mean void size does not 
change much. We get a value of D mean ~ 10 Mpc (18 Mpc) in a 64 3 sample and 
D mean ~ 10.5 Mpc (19 Mpc) in a 128 3 sample. (The numbers quoted in paren- 
theses are those obtained after discarding all voids of size less than 10 Mpc as in 
KF.) This could indicate that the mean diameter of voids converges rapidly in an 
= 1 CDM Universe so that the results which we quote for D mean from our 128 3 
simulations could be representative of the entire CDM Universe. However, this is- 
sue can be completely resolved only after larger simulations have been performed. 
(Incidentally, the maximum void size increases from -D max — 40 ± 8 Mpc in a 64 3 
simulation to -D max — 57 ± 6 Mpc in a 128 3 simulation.) We would also like to 
mention in this connection, that the median (evaluated with respect to the void 
number spectrum) is usually a considerably more stable statistical indicator of void 
sizes than the mean. The reason for this is that the median is much less sensitive to 
the presence of a few very large voids in a sample than the mean, and consequently 
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its value does not grow as rapidly as the latter as we increase the sample size of our 
model Universe. (See also Kauffmann & Melott 1992). 

Fig. 1 shows that one might expect to see a correlation between the sizes of 
voids, and the height of the linear gravitational potential, especially at late times. 
This feeling is borne out by the results of our simulations for the CDM model (see 
Fig. 5), which show that a distinct correlation exists between the size of a void, 
and the value of the gravitational potential evaluated at its center. The value of 
the potential averaged over all the void centers at the present epoch is about 0.75 
times the rms value of the potential, indicating that voids preferentially form in 
those regions where the gravitational potential is large. In order to demonstrate 
the general nature of this relationship between void size and the value of the grav- 
itational potential 0, we have also analyzed simulations for power law primordial 
spectra: P(k) = (\5k\ ) oc /c n , where \Sk\ is the spatial Fourier transform of the 
density contrast (see Melott & Shandarin 1993). It is assumed that the phases of 
the Fourier components are randomly distributed, so that the gravitational poten- 
tial (f)(x) has the statistical properties of a Gaussian random field. A cutoff was 
introduced into the spectrum by requiring that P{k) = for k > k c . (Physical 
processes which can give rise to such a cutoff include the free streaming of weakly 
interacting massive particles (WIMP's) such as massive neutrino's, and the Silk 
damping of inhomogeneities in the photon-baryon plasma.) Our simulations were 
run on a 128 3 box for the following three models: n = —1, 0, +1, and k c = 16 x kf. 
(kf is the wavenumber of the fundamental mode: kf = 2it/L where L is the length 
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of the side of the simulation cube.) Three simulations were performed for each of 
the three spectra amounting to 9 simulations in all. 

Our results are shown in Fig. 6 for different values of cr(t). The upper limit in 
the integral in equation (10) is effectively replaced by k^ - the Nyquist frequency 
in the simulation. a(t) can be related to k^L ~ the wavenumber of a mode going 
nonlinear at the given epoch: 

»p/p) 2 ) = ^y o p(k)k 2 dk = i (12) 

so that 

a(t) = ^ . 13 

\ft»<-P{k)k*dkJ 

We see that for truncated power law spectra (n > —3) 

consequently, k^L decreases with epoch, as successively larger scales enter the non- 
linear regime. We do not evolve our simulations beyond the epoch when the largest 
scale to go nonlinear equals the box size. 

From Fig. 6 we find that as in the case of the CDM model, a marked correlation 
exists between the diameter of a void and the value of the gravitational potential at 
its centre at late times. We also find that the mean value of the potential evaluated 
at void centers 

1 N 
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N being the number of voids in a simulation, increases with time, being larger at a 
given instant of time, for steeper spectra (see Fig. 7). (4>)/4> r ms oc y/cr(t) provides 
a reasonably good approximation, especially at late times, to the rate of growth of 
((f)) for the power law spectra considered by us. 

From figures 5, 6 we find that the correlation between diameter and potential is 
more striking for spectra with more short range power (such as the truncated n = 1 
spectrum). This may be because a potential having significant long range power 
(such as CDM) has features on it (small scale wiggles superimposed on long range 
mountains and valleys) which will influence the sites of void formation. As a result 
a small wiggle in the potential will generally give rise to a small void even if the 
wiggle occurred at a large value of ((f)) / '(f) rms . (The dip in ((f)) /(f) rms occuring for void 
dimaters 10 — 20 Mpc. in figures 5, 6, which is more pronounced for spectra with 
more large scale power, is a distinctive feature of this effect.) This modulation in the 
void spectrum caused by low frequency waves in the potential, will lead to a more 
noisy relation between the void diamter and (4>)/4>rms for spectra with significant 
large scale power, as is indeed demonstrated in figures 5,6. 

The correlation of void sizes with the "height" of the gravitational potential 
seen in Fig. 5 and 6, raises the interesting possibility of reconstructing the primordial 
form of (f) from the observed form of the void spectrum. The primordial form of 
the gravitational potential is determined both by physical processes occuring in the 
very early Universe (such as Inflation), as well as by the nature of dark matter, 
hence some indications as to its form would be of immense value (see also Kofman 
& Shandarin 1988). 
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The evolution of the void spectrum and the transfer of power from smaller to 
larger voids at late times which was shown to occur for the CDM model is demon- 
strated with great clarity in Fig. 8 for power law spectra. (Note particularly the 
evolution of the void number histogram shown in Fig. 8 b-d and the corresponding 
evolution of the mean and maximum void sizes in Fig. 9a.) We find that for large 
voids V > Vmean the number spectrum of voids in Fig. 8b - 8d is well fitted by the 
exponentially decaying function exp(— V/V*) where V* oc V mean . A similar result 
was also obtained by Kofman et al. (1992) for the two dimensional case. 

From Fig. 9a we find that the mean void diameter grows very nearly as D mean oc 
a/ a(t) at late times. (The complementary picture demonstrating the evolution of 
the total number of voids in our simulations is shown in Fig. 9b.) For the truncated 
scale invariant spectrum n = 1, this result agrees with an asymptotic analysis of 
the growth of cellular structure in the adhesion model made by Gurbatov et al. 
(1985, 1989), who found that D mean oc yfifi) if n > 1, and D mean oc a(t) 2 /( n+Ar °) 
(Nd being the dimensionality of space) for — 1 < n < 1. The last result suggests 
that the mean void size grows faster for spectra with more large scale power, and 
is seemingly at odds with our analysis which indicates that D mean oc a/ a(t) for 
— 1 < n < 1. This discrepency might arise because of the following two reasons: 
a) In their analysis Gurbatov et al. implicitly assume that voids are associated 
with peaks in the gravitational potential. We find on the other hand, that this 
assumption is not always satisfied, since voids can sometimes be associated with 
other features in the potential such as ridges etc. (this is especially true in two 
and three dimensions.) As a result the analysis of Gurbatov et al. is strictly valid 
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at asymptotically late times, b) Our simulations are carried out in a finite sized 
box which might inhibit the late time growth of voids particularly for fluctuation 
spectra with long range power. A similar effect was earlier discussed by Kofman et 
al. (f 992) in the context of the two dimensional adhesion model. Indeed, simulations 
conducted by us on spectra having no small scale cutoff show that D mean grows 
proportionally to a(t) 2 /( n+ND "> for n = 0, — 1 as expected. 

An interesting feature which emerges from our simulations with a cutoff in the 
initial spectrum is that the number of voids peaks at a given value of the expansion 
factor a(t) = a*, and then declines steadily (see Fig. 9b). The value of a* is 
relatively insensitive to the form of the spectrum for spectra with k c ^> kf, and can 
be described by the analytic formula: 

C* = 0"(&7VL = R*) (15) 

where R* = v / 3cn/o"2- &j are the moments of the distribution of <fi (Bardeen et al. 
1986, Kofman 1989): 

°i = 2^2/ dkk2k2j \^\ 2 - (16) 

The scale R* characterises the mean distance between peaks of the gravitational 
potential. For the truncated power law spectra considered by us: k c x R* ~ 
\/6, 3, and 4, for n = 1, 0, and -1, respectively. For spectra having no built-in 
scale, the value of R* would be determined by the small scale cutoff introduced by 
the Nyquist frequency in a simulation.) 

Fig. 10 shows the fraction of matter in caustics (stuck matter) f(a) (dotted 
curves), and the related underdensity of matter in voids 1 — f(a) (solid curves) as 
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a function of epoch, for the same set of simulations. Our simulations show that the 
rate of infall of matter into caustics is virtually insensitive to the spectral index in 
the range — 1 < n < +1. We also find that by a(t) = a*, virtually all of the matter 
is already in caustics, so that voids are relatively empty by this epoch (see also 
Pogosyan 1989, Sahni 1991). (The fraction of matter in caustics can be described 
to within 2% accuracy by the fitting function f(a) = a 3 6 (3 + a 3 ' 6 ) -1 for a > 1.) 

Fig. 9b lends support to the viewpoint that the growth of large scale structure 
in the Universe is characterised by two complementary regimes. The first regime 
describes the formation of pancakes and the establishment of cellular structure 
(cr(t) < o* during this period). However, the emergence of cellular structure in 
the Universe, is only an intermediate asymptote, and is succeeded by a second 
epoch during which large mass concentrations (knots and filaments) attract one 
another leading to the disruption of cellular structure at late times. As a result, 
whereas gravitational instability leads to pancaking during ait) < a*, hierarchical 
clustering characterises the later period a(t) > a*. This interpretation is confirmed 
by a comparison of the results of 2D N-body simulations (Fig. 11 a-c) (Melott & 
Shandarin 1989; Beacom et al. 1991) with those of the adhesion model (Fig. 11 
d-f). From Fig. 11 it is clear that the formation of cellular structure is complete 
by a = a* ~ 4 and that from then onwards clustering proceeds hierarchically. In 
fact, from Fig. 11 a-c and d-f we find that the number of voids first increases from 
one interconnecting void at a = 1, to a maximum of eight voids at a = 4 and then 
decreases to five voids at a = 16. (We would like to add that despite the fact that 
gravitational clustering proceeds hierarchically during the second epoch (a(t) > cr*), 
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we can still find the emergence of large correlated structures such as secondary 
pancakes, during this period, as has recently been demonstrated by Kofman et al. 
1992.) This view of the formation and evolution of large scale structure is borne out 
by a study of the evolution of the void spectrum shown in Fig. 4 and Fig. 8a-d. From 
these figures we find that there are a greater number of large voids at early and late 
times. This clearly demonstrates that the evolution of voids proceeds in two distinct 
phases. During the first phase, (u(t) < a*) large voids fragment into smaller voids. 
During the second phase, (a(t) > a*), the reverse process takes place as voids begin 
to coalesce into larger units. The mean and maximum void sizes shown as functions 
of a(t) in Fig. 9a, provide further support to this point of view. Strictly speaking 
the above model assumes a cutoff in the initial spectrum as in our simulations. 
However, if the initial spectrum falls down as /c~ 3 , as in the CDM model, then 
the evolution of the large scale structure also follows the above discription. The 
first caustics in this case are associated with the free-streaming distance for dark 
matter particles - in the case of a real (spatially unbounded) Universe, or with the 
artificial, Nyquist frequency cutoff - for a numerical simulation. 

Finally, we would like to comment on the fact that the correlation of the height 
of the primordial potential with the void size appears to be more noticeable at late 
times when a > cr* (see Fig. 5 and 6). The reason again lies in the fact that the 
formation of cellular structure is complete only by a ~ cr*, where a* ~ 5 for the 
power law spectra discussed in the previous paragraph. Therefore at earlier times 
when a < cr* large voids are still in a state of becoming, since what used to be a large 
void at an early epoch gets fragmented into several smaller voids at later times and 
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vice- versa. As a result, at epochs prior to a ~ a*, no definite relation as yet exists 
between the sizes of large but transient voids, and the height of the gravitational 
potential within them. At very late times the gravitational potential at the centres 
of voids is very nearly uniform irrespective of their size. (We should also note 
that since there are very few large voids in a given simulation, their distribution is 
plagued by small number statistics as evident from Fig. 5 and 6.) 



A central result of our study of voids is that voids can be populated by sub- 
structures such as minipancakes and filaments which run through a void. Tis result 
emerges from a study of the topology of voids in Lagrangian space. The topology of 
a compact manifold, can be characterised by its genus measure g which is related to 
the integrated Gaussian curvature K of the manifold, by the Gauss-Bonnet formula 



(For a sphere K = r~ 2 ,A = 4nr 2 , and therefore g = 0.) The genus is related to 
the Euler characteristic of a manifold k, by the simple relation g = 1 — k/2. At 
an intuitive level, we might say that the genus characteristic provides a measure 
of the number of holes in a given manifold. Thus, a sphere has genus = 0, a 
torus genus = 1, a pretzel genus = 2, and so on. In our present discussion we 
shall be interested in evaluating the topology of voids defined on a discrete grid in 



IV. SUBSTRUCTURE WITHIN VOIDS 




(17a) 



30 



Lagrangian space. To do this we shall find it convenient to work with the discrete 
analog of the Gauss-Bonnet formula: 



Here Di is the angle deficit defined at each vertex of the polyhedral void. Equation 
(17b) can be derived from equation (17a) if we note that all of the curvature of a 
discrete surface is localised in each of its N vertices so that 



We apply equation (17b) to determine the topology of individual voids defined in 
L-space using the adhesion ansatz. (We would like to stress that, throughout the 
present discussion, the topology that is referred to will be that of individual voids, 
as opposed to the topology characterising large scale structure as a whole. See 
Melott 1990 for a comprehensive review of the subject of the topology of large scale 
structure.) 

One of the most important results to emerge from the present analysis of voids 
using the adhesion model, is that voids in L-space can have a nontrivial topol- 
ogy. (We would like to point out in this context that the 2 dimensional L-space 
voids shown in Fig. 2c and Fig. 11 e,f, are all homomorphic to circles and have 
therefore an essentially trivial S 1 topology. In contrast to this the voids appear- 
ing in our three dimensional simulations sometimes show appreciable departures 




(176) 



K{x) ~ ^Di8{x-Xi) 



and 
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from a purely spherical (i.e., S 2 ) topology.) To illustrate what we mean, we plot 
one of the several topologically nontrivial void configurations that appear in our 
simulations in Fig. 12. The void has the topology of a torus in L-space i.e., the 
free region defining the void has a single hole (corresponding to stuck particles) 
passing through it. As mentioned in section 2 of this paper, the Euler picture of 
the distribution of caustics can be constructed from the coresponding Lagrangian 
picture, by moving the border between stuck and free regions by means of the Zel- 
dovich approximation. This essentially amounts to moving the boundary of a void 
in L-space by the Zeldovich prescription (1). Since the Zeldovich approximation 
is a topology preserving transformation (Shandarin & Zeldovich 1989), it follows 
that a non-simply connected void in L-space will be mapped onto a topologically 
non-trivial configuration in E-space. In particular, a void having the topology of 
a torus (i.e., genus = 1), will correspond to a void having one minor Zeldovich 
pancake running through it. In general, a void in L-space having genus = N g , will 
correspond in E-space, to N g minor Zeldovich pancakes running through a void 
bounded by major Zeldovich pancakes. (Observationally, such voids might convey 
the impression of being sponge-like.) We also feel that mini-pancakes within voids, 
could be young transient features, since our simulations show that voids tend to 
empty out with time, leading to a bubble like topology for voids at late times. The 
dynamical mechanism leading to the trivialisation of the void topology could be the 
following: As voids expand they empty out with time and matter from within them 
falls into pancakes (both mini as well as major). As a result, mini pancakes within 
voids evolve with time growing progressively more massive. At a later epoch these 
minor pancakes graduate into major Zeldovich pancakes and succeed in dynamically 
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dividing a topologically nontrivial void into two or more topologically trivial voids. 
In order to estimate the evolution of the void topology with time, we determine the 
mean genus characteristic of all voids V\ in a given simulation 



N being the total number of voids at a given epoch. 

The results of our simulations for the CDM model, are shown in Fig. 13 and 
14. In Fig. 13 we plot the histogram of the genus per void against void size for three 
different epochs. What is plotted is (gi) - the total genus of all voids in a certain 
class divided by the number of voids in that class: 



where is the number of voids in the diameter class Di and g^, k = 1, . . . , are 
the genus values of the individual voids in that class. This is a good estimator of 
genus as a function of void diameter since the number of voids in different classes are 
different. In Fig. 14 we plot the mean genus in our simulation as a function of <j(t). 
These two figures reveal that one out of every fourteen voids today has a nontrivial 
topology. By contrast, every eighth void had a topology that was nontrivial at 
a redshift of unity. Fig. 14 clearly demonstrates how the void topology evolves 
and essentially trivialises with time. A good analytical approximation to the time 
evolution of the mean genus, shown in Fig. 14 for the CDM model, is provided by 
g(a) oc a~ 3 / 2 (t). For comparison we also plot the genus, against void size and as a 
function of cr(t), for the power law spectra P(k) = (\5k\ ) oc k n for k < k c , P{k) = 




(18) 



i=l 
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for k > k c , n = — 1,0, +1, in Fig. 15, & 16, respectively. We find that in this 
case too, larger voids tend to have a more complicated topology than smaller voids, 
and that the overall void topology tends to simplify with time. It is interesting to 
note that some voids such as the Coma void, do indeed seem to be populated by 
mini-Zeldovich pancakes, as pointed out recently by Park et al. 1992. 

From Fig. 13 we find that the ratio of the diameter of a void in Eulerian space 
to its diameter in Lagrangian space, grows very nearly as De/Dl oc cr(t), regardless 
of the size of the void. This relation is easy to understand if one notes that at late 
times equation (11) reduces to 

N 

V E =J2dV L [1 - a(t)A!( gi )] [1 - a(t)X 2 ( qi )} [1 - a(t)X 3 ( qi )\ 

(19) 

= -Ndv L a\t) y hMhMhM 

a»l ^— ' N 

so that, 

V-J =^ = *(')</s>', 
1 N 

(h) = -^Yl A i(9<) A 2(ft)A 3 (gi), 
i=i 

is the mean value of the invariant 7 3 in the given L-space void. iV denotes the 
number of free unit volume elements in a given L-space void, and Vl = NxdV^ is the 
Lagrangian volume of the void (cr(t) oc a(t)). From Fig. 13 it follows that small voids 
in L-space remain small in E-space and vice versa. The ratio Vl/Ve = (D L /D E ) 3 
is a measure of the underdensity of voids. From this figure we find that 

Vl\ _\/Yl\ 
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where 



indicating that voids today are ~ 8 times more empty than they were at redshifts 
of unity. (Note that a ~ 9 today, for the COBE-normalised CDM model.) 

V. DISCUSSION 

We have applied the adhesion model to determine the spectrum of void sizes, 
both for the COBE-normalised CDM model, as well as for models with power-law 
spectra P(k) oc k n , n = — 1,0, +1. The model neglects the volume occupied by 
clusters and superclusters, however it does not make any a priori assumptions on 
the geometry and topology of the voids. The advantage of the adhesion model 
is that it includes both dynamical as well as statistical aspects of void formation. 
We find that most of the characteristics of voids, such as the mean and maximum 
void size in a simulation, evolve with time. For the CDM model we find that 
D = Do/y/1 + z, where Do — 10.5 Mpc, is the mean void size today. We also find 
that the total number of voids in a simulation, varies with cosmic epoch, increasing 
to a maximum value at a* = a(k^ f L L = R*), and then decreasing monotonically. 
We feel that the expansion scale cr* separates two distinct phases in gravitational 
instability. During the first phase (cr < a*), gravitational instability proceeds via 
the formation of pancakes, and the subsequent establishment of cellular structure, 
which is complete by a(t) = a*. However, pancake formation happens to be an 
intermediate asymptote and during the second phase (a(t) > a*) gravitational 
instability proceeds hierarchically, leading to the eventual disruption of cellular 
structure. 
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Our results show that the sizes of voids are stongly correlated with the height 
of the primordial gravitational potential at void centers, larger voids forming in re- 
gions where the gravitational potential is higher. This result indicates that knowing 
the observed form of the void spectrum it might be possible to, at least partially, 
reconstruct the primordial form of the potential. Assuming the Harrison- Zeldovich 
initial spectrum one can easily verify that the gravitational potential is not a homo- 
geneous random function if the spectrum is extrapolated to k = 0. It means that 
the largest void is probably determined by the size of a sample. 

One of the most intriguing results of the present analysis is that voids can have 
a topology which is nontrivial. We find that the topology of a void depends upon 
its size, with larger voids more likely to have a nontrivial topology than smaller 
voids. We also find that the topology of voids, like the void spectrum evolves with 
time, with voids becoming progressively emptier at later times. Our results for the 
CDM model show that one out of every 14 voids is likely to have some substructure, 
such as a filament or a pancake passing through it. These results are supported by 
recent observations indicating the presence of galaxies (Westrop et al. 1992) and 
mini-pancakes (Park et al. 1992) within at least two voids. It is likely that with 
the advent of deeper and more complete redshift surveys the issue of void size and 
void substructure will take on a deeper significance, as a detailed picture of the 
texture of voids in the Universe emerges. Bearing this in mind, we believe that the 
topological indicator of void substructure developed by us in this paper, can emerge 
as a key statistical indicator for the study of voids in the Universe. 
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If both the shape and the amplitude of the initial perturbation spectrum are 
independantly known (for instance from the angular fluctuations of the cosmic mi- 
crowave background radiation), then void statistics in the adhesion approximation 
is determined solely by the growth factor of the linear density contrast D(z,Q) 
(which determines the curvature of the osculating paraboloid). In such a case the 
void spectrum can potentially be used to determine the value of O. 

In a companion paper we shall apply the void spectrum and topology measuring 
ansatz, developed in the present paper, to compare and contrast several cosmological 
scenario's, including "tilted" cold dark matter models, cold dark matter models with 
a cosmological constant, and models containing a mixture of cold + hot dark matter 
(Sahni, Sathyaprakash & Shandarin 1993). 

Finally we would like to mention that the consistency of the results obtained 
in this paper was checked using the 3D adhesion code independently developed by 
Dima Pogosyan. 
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FIGURE CAPTIONS 



Fig. 1 

The geometrical prescription of descending a paraboloid onto the gravitational 
potential in order to demarcate stuck and free Lagrangian regions, is shown in one 
dimension. The peaks of the potential correspond to regions where the eigenvalue 
A < 0. The particle having Lagrangian coordinate qo is free in the uppermost 
figure, and has just entered into a caustic in the middle figure. The middle and 
lower figures describe the merger of clumps. 

Fig. 2 

The distribution of caustics (dots) (for two-dimensional simulations) is plotted 
superimposed on free (shaded) and stuck (unshaded) Lagrangian regions for three 
distinct expansion epochs: (a) a(t) = 1, (b) a(t) = 2, and (c) a = 8. (d) Results 
of two-dimensional N-body simulations are shown for the epoch a(t) = 2 (Melott 
& Shandarin 1989). (e) Particles evolved according to the Zeldovich approximation 
are shown for a(t) = 2. 

Fig. 3 

Regions in Lagrangian space where both the eigenvalues of the two-dimensional 
deformation tensor are negative, are shown (dotted regions). Also shown superim- 
posed are voids in Lagrangian space obtained using the adhesion model, for the 
same spectrum, at very late times (shaded regions). 
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Fig. 4 



The void spectrum for the COBE - normalised CDM model, is plotted for three 
time scales: (top) the past a = 0.5 ao, z = 1; (center) the present a = ao, z = 0; and 
(bottom) the future (!) a = 1.5 ao; (a - is the expansion factor, and z the redshift). 
In all the pictures the x-axis is an indicator of the void diameter in Mpc. The left 
hand pictures illustrate the volume fraction of voids having a given diameter "£)" 
(in Mpc) v(D)/V sometimes called the void spectrum. (v(D) is the volume occupied 
by voids in a bin of size D ± 2 Mpc and V is the total volume occupied by all voids 
which, in the adhesion model, is equal to the volume of the simulation box.) The 
right hand pictures show the corresponding number fraction of voids, also plotted 
against the void diameter. The error bars in these histograms correspond to rms 
dispersion over three simulations. We would like to point out that the number of 
large voids in a given simulation is too small to be resolved in the number fraction 
histogram. Their presence is, however, seen in the void spectrum. 

Fig. 5 

The value of the normalised primordial gravitational potential (0/0 rms ) eval- 
uated at the centers of voids is shown plotted against the void diameter for the 
present epoch for one realisation of the COBE - normalised CDM spectrum. 

Fig. 6 

The value of the normalised primordial gravitational potential (4>/4> r ms) eval- 
uated at the centers of voids is shown plotted against the void diameter for power 
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law primordial spectra P(k) oc /c n , for k < k c ; P(k) = 0, for k > k c , where 
k c = 16 x kf, kf = 2ir/Lbox, is the wavenumber of the fundamental mode. The fig- 
ures (top to bottom) correspond to three different cosmic expansion epochs: a = 4, 
a = 8, and a = 16. From left to right the figures correspond to n = 1, (left column) 
n = (middle) and n = — 1 (right). 

Fig. 7 

The growth of ((f)(t)) / 4> rms is shown plotted against the cosmic expansion factor 
cr(t), for two of the three power law spectra discussed in the previous figure: n = 1 
(solid line) and n = (dashed line). 

Fig. 8 

(a) The void spectrum v(D)/V is plotted against the void diameter D, for the 
truncated power law spectrum of Fig. 6, 7 with spectral index n = 1 and for nine 
different expansion epochs, (b) The number fraction of voids n(D)/N is plotted 
against the void diameter D for the power law spectrum and expansion epochs 
of Fig. 8a. (c) The void spectrum v(D)/V (left) and the void number fraction 
n(D)/N (right) are plotted against the void diameter D for the truncated power 
law spectrum of Fig. 6, 7 with spectral index n = and for three different expansion 
epochs, (d) The same as Fig. 8c but with spectral index n = — 1. 

Fig. 9 
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(a) The time evolution of the mean (left panel) and maximum (right panel) 
void diameters, in our simulation cube, are shown as functions of a(t) for the same 
power law primordial spectra as in the previous figure: n = 1 (solid line), n = 
(dashed line), n = — 1 (dotted line). The error bars indicate the rms dispersion 
from three simulations, (b) The time evolution of the total number of voids in our 
simulation cube, is shown as a function of a(t) for the same power law primordial 
spectra as in the previous figure: n = 1 (solid line), n = (dashed line), n = — 1 
(dotted line). The error bars indicate the rms dispersion from three simulations. 

Fig. 10 

The fraction of matter in caustics (dotted curves) and the corresponding un- 
derdensity of matter in voids (solid curves) are shown as functions of the expansion 
factor cr(t), for the power law primordial spectra previously considered. (The up- 
per (lower) dotted curve and the lower (upper) solid curve correspond to n = — 1 
(n = 1).) The dispersion resulting from different simulations is so small that it 
cannot be resolved in this graph. Also, the curves are virtually insensitive to the 
index of the power spectrum. 

Fig. 11 

Results for two-dimensional N-body simulations involving 512 x 512 particles, 
are shown for different instants of time corresponding to: (a) a(t) = 1, (b) a(t) = 4 
and (c) a(t) = 16. The primordial power spectrum for these simulations was 
assumed to be a truncated power law: P(k) oc k 2 for k < k c , P(k) = for k > k c , 
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where k c = 4 x kf, kf being the fundamental frequency kf = 2ti/L} }OX . (Melott & 
Shandarin 1989.) Voids in L-space (described by shaded regions) evaluated using 
the adhesion approximation are shown for initial conditions and epochs identical to 
those in Fig. 11 a-c: (d) a(t) = 1, (e) a(t) = 4, and (f) a(t) = 16. 

Fig. 12 

A topologically nontrivial void in Lagrangian space having genus = 1 is shown 
from our CDM simulations. 

Fig. 13 

In the left hand pictures, the mean genus of voids belonging to a given diameter 
class (gi), is shown plotted against the void diameter in Lagrangian space Dl, 
for simulations involving the CDM model, and for the same expansion epochs as 
in Figure 5 (i.e., from top to bottom: a/ao = 0.5, 1, 1.5). The corresponding Euler 
diameters of voids, De, are shown plotted against the Lagrangian diameters Dl, in 
the right hand pictures. The solid line corresponds to De/Dl = cr{t)(\\i\2^3\)^ , 
which is a constant for a given simulation. 

Fig. 14 

The mean genus characteristic evaluated for the entire ensamble of voids in the 
CDM model, is shown for the expansion epochs: a/ao = 0.5, 1, 1.5, corresponding 
to the past {z = 1), the present (z = 0), and the future. 

Fig. 15 
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Same as Fig. 13 but for power law spectra P(k) oc k n , k < k c ; P(k) = 0, k > k c , 
where k c = 16 x kf, kf = 27r/Lb ox , is the wavenumber of the fundamental mode: 
(a) n = 1, (b) n = 0, and (c) n = —1. 

Fig. 16 

Same as Fig. 14 but for the power law spectra of Fig. 15. 
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